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ABSTRACT 

The start— up process of a frozen heat pipe is 
described and a complete mathematical model for the 
start-up of the frozen heat pipe is developed based on the 
existing experimental data, which is simplified and solved 
numerically. The two-dimensional transient model for the 
wall and wick is coupled with the one-dimensional transient 
model for the vapor flow when vaporization and 
condensation occur at the interface. A parametric study is 
performed to examine the effect of the boundary 
specification at the surface of the outer wall on the 
successful start— up from the frozen state. For successful 
start-up, the boundary specification at the outer wall 
surface must melt the working substance in the condenser 
before ary-out takes place in the evaporator, 

NOMENCLATURE 

Cp specific heat at constant pressure, J /kg— K 
c v specific heat at constant volume, J/kg— K 

D diameter of the vapor space, m 
£ total energy of the vapor per unit volume, p [c v T y 4 * 

1/2 (V 2 y + Vy)j 

f friction coefficient at the wall, r/pU^ 

h enthalpy, J/kg 

h c heat transfer coefficient, W/m - K 
h^ latent heat of vaporization, J/kg 

o 

latent heat of melting per unit volume, J/m 

2 

k permeability of the wick structure, m 

K thermal conductivity, W/m-K 
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L length of the heat pipe, m 

m Q rate of evaporation or condensation per unit area, 

2 

kg/m -s 

M molecular weight, kg/Kmol 

n unit outward normal direction 

o 

P pressure, N/m~ 

2 

q heat flux, W/m 

9 

q' new heat dux, W/m“ 

2 

q" new' guessed heat finx, W/m 

P cr reference pressure for the Clausius-Clapeyron 
*) 

relationship, N/m" 

r coordinate in the radial direction 

r c effective capillary radius, m 

ip radius of curvature of the meniscus, m 

Ar radial distance between nodes, m 
R q universal gas constant, KJ/Kmol— K 

R y radius of the vapor space, m 

R w inner radius of the heat pipe wall, m 

S interface position in vector, m 

t time, s 

T temperature, K 

T c reference temperature for convection, K 
T^ reference temperature for the Clausius-Clapeyron 
relationship, K 

T m melting temperature, K 

T r reference temperature for radiation, K 

T s saturation temperature, K 

Tj_^ temperature at node j-1 in the radial direction in the 
wick region, K 

Tj temperature at last node near liquid— vapor interface 
^ in the radial direction in the wick region, K 
T transition temperature, K 
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AT small finite temperature interval around T to 

m 

define mushy zone, K 
U axial velocity, m/s 

V radial velocity, m/s 

x coordinate in the axial direction 

GREEK SYMBOLS 
a relaxation factor 

4(T~T m ) Dirac function 

c porosity 

8 contact angle of the liquid, deg 
A length of mean free path, m 

o 

p dynamic viscosity, N— s/m 

2 

u kinematic viscosity, m /s 

f emissivity 

p density, kg/m 

0 A 

a Stefan-Bohzmaan constant, W/(nT— K ) 

r shear stress 

ij surface tension, N/m 

SUBSCRIPTS 

it working substance in the liquid state in the wick 

fs working substance in the solid state in the wick 

/ liquid where there is liquid motion in the wick 
te wick region where the working substance is in the 
liquid state 

me wick region where the working substance is in the 
mushy state 

o properties at the liquid— vapor interface 

s wick structure material 

se wick region where the working substance is in the 
solid state 
v vapor state 

w properties at the heat pipe wall 

INTRODUCTION 

The demand for an effective thermal management 
device for high temperature applications such as cooling the 
leading edges of reentry vehicles and hypersonic aircraft, 
and a space-based power station stimulates the study of 
the start— up of frozen liquid— metal heat pipes. Also/ the 
start-up of frozen low temperature heat pipes is important 
in applications such as heat pipe heat exchangers, cooling 
electronic equipment, and melting the snow and ice on roads 
and bridges. 

Neal (1967) and Shiosinger (1968) performed the first 
experimental tests to study the start— up performance of low 
temperature heat pipes with the water initially frozen. 
Deverall et al. (1970) also made a series of tests with wafer 
and liquid metal heat pipes. Successful start— up from the 
frozen state was possible but was highly dependent on the 
heat rejection rate at the condenser. Tolubinsky et al. 
(1978) investigated the start-up characteristics of sodium 
and potassium heat pipes. Camarda (1977) investigated the 
performance of a sodium heat pipe cooling a leading edge. 
Start— up and shut-down of a 4 m long lithium heat pipe 
was studied experimentally by Merrigan et al. (1985, 1986). 
Ivanovskii et al. (1982) presented the vapor temperature 
distribution along the length of a sodium heat pipe during 
the start-up period. Three flow regimes in the condenser 
are described based on the vapor temperature: free 
molecular flow, intermediate, and continuum vapor flow. 
Unfortunately, the existing experimental data for the 
start-up period are in general not represented in the 
archival literature and lack sufficient information for 
comparison with numerical simulations. 

Colwell et al. (1987) and Jang (1988) developed a 
simple mathematical model to predict the start-up behavior 
of a sodium heat pipe with a rectangular cross section from 


the frozen state. In the wall and wick structure. energy 
transport is described by the transient, two-^nmeasional 
heat conduction equation, and the phase change of the' 
working substance is taken into account. In the vapor 
region, free molecular, choked and continuum flows are 
considered and one— dimensional, compressible quasi-steady 
state laminar flow is assumed. The numerical results 
obtained by using the finite element method are in 
agreement with experimental results given by Camarda 
(1977). 

The liquid metal heat pipe operates not only at high 
temperatures but also the initial temperature may be 
ambient temperature. In this temperature range, the 
working substance may be in the solid state as well as the 
liquid and vapor states. In the vapor space, free molecular 
flow, continuum flow, sonic and supersonic flow may be 
encountered due to the extremely small density during the 
start-up of the heat pipe. These conditions may cause the 
failure of operation of the heat pipe and limit the 
performance of the heat pipe. Understanding the start— up 
behavior and transient performance of the high temperature 
heat pipe is therefore important and an efficient 
mathematical model is needed to predict this behavior. 

The first part of this paper presents a complete 
mathematical model to describe the start-up behavior of 
the heat pipe from the frozen state. This model is then 
simplified to obtain numerical results. To the authors 5 
knowledge, the analysis presented here is the oniy model 
that includes the effect of the transient vapor flow in the 
analysis of the start-up of a heat pipe from the frozen state. 
After the mathematical models are tested separately, the 
models are used for a parametric study of the start-up of 
frozen heat pipes. 

DESCRIPTION OF HEAT PIPE START-UP 

Previous experimental observations suggest the 
following sequence of events during heat pipe start— up from 
the frozen state. Initially, the working substance is in the 
solid state and the vapor density is extremely low, so that 
free molecular flow conditions prevail throughout the vapor 
space. The input heat flux over the evaporator increases the 
temperature and starts to melt the frozen substance in this 
region. Meanwhile, the heat transport from the heated zone 
to the adjacent pipe proceeds quite slowly via axial 
conduction through the heat pipe wall, working substance, 
and wick structure, while the heat transfer in the vapor is 
almost negligible. Thus, a large temperature gradient exists 
between the evaporator and condenser. 

When energy is continuously added to the 
evaporator, the frozen working substance in the evaporator 
is melted, so that evaporation can take place at the 
liquid— vapor interface and the vapor density in this region 
is increased. The molecular mean free path in the heated 
region then becomes small compared to the diameter of the 
vapor passage and the continuum flow' regime is established, 
while in the cooled zone the vapor is still in free molecular 
flow. In the continuum flow region, the vapor flows into the 
condenser section due to the large pressure gradient. During 
this stage, energy is mainly transferred as latent heat owing 
to vaporization in the heated zone, and condensation in the 
cooled zone in the vapor space where continuum flow is 
established. The temperature near the evaporator remains 
constant and the location of the temperature gradient moves 
toward the end of the condenser until continuum flow is 
establised in the entire vapor space. Cotter (1967) also 
described this frontal start-up mode when the vapor density 
is so low and non-con densable gas does not exist in the 
vapor space. 

In heat pipes with metallic working substances, the 
vapor densities are very small during the start-up even in 
the continuum flow region. Thus, even for relatively low 
values of heat input, sonic vapor velocities can be reached. 
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Also, the vapor flow in the heat pipe is quite similar to the 
flow in a conver^ing^di verging nozzle due to the vapor 
addition in the evaporator and the vapor removal in the 
condenser (Dunn and Reay, 1982). Thus, the heat transfer 
through the vapor space may be limited by the choked flow- 
condition, and supersonic vapor flow and a shock front may 
occur in the continuum flow region in the condenser. The 
maximum rate of heat transfer is limited by the sonic limit 
so that a high heat input in the evaporator causes the 
various types of start-up failure. 

The start-up process continues until the frozen 
working substance is completely melted and the continuum 
flow regime reaches the end of the heat pipe, at which time 
liquid returning to the evaporator is sufficient for normal 
transient operation. Eventually, the heat pipe may reach a 
steady state condition. The start— up process of the liquid 
metal heat pipe from a frozen state may be divided into 
several distinct periods for convenience of analysis based on 
the status of the working substance and the behavior of the 
vapor flow. 

1. In the first period, no phase ch*nge takes place in the 
entire region but the temperature near the heated 
region increases. The vapor flow is in the tee 
molecular condition. 

2. The working substance in the evaporator is in the 
liquid state, bur evaporation does not occur at the 
liquid— vapor interface. 

3. The liquid and solid states of the working substance 
exist simultaneously in the wick structure and 
vaporization of the working substance takes place at 
the liquid— vapor interface. In the vapor space, a 
region of continuum flow is established in the heated 
region and a continuum flow front moves iow-ard the 
cooled end of the heat pipe. The vapor flow may be 
choked at the beginning ox the condenser. 

4. The working substance is completely melted but free 
molecular flow still exists in part of "the vapor space. 

5. Continuum flow exists over the entire heat pipe 
length in the vapor region but the heat pipe does not 
reach the steady state condition. 

5. The heat pipe then reaches the steady state 

operation. 

For low temperature heat pipes, the experimental 
resuits of the successful start— up from the frozen state are 
very rare. Deverail et ai. (1970) successfully started a water 
heat pipe from the frozen state (208 K). The wall 
temperature distribution obtained is similar to that of high 
temperature heat pipes. The vapor temperature was not 
obtained, but the vapor density is relatively high even 
around the melting temperature. This means that the vapor 
velocity is very low so that choked flow and supersonic 
vapor velocities may not be encountered during start-up. 
Unlike high temperature heat pipes, experimental results 
show that the heat pipe becomes immediately active where 
the ice is melted, but there still is a large temperature 
gradient in the axial direction. The heat input at the 
evaporator should be low enough to return sufficient water 
into the evaporator for the successful start-up. 

MATHEMATICAL FORMULATION 

Consideration is given to the heat pipe wall, the wick 
structure with the working substance initially subcooled at a 
uniform temperature, and the vapoT region. A schematic 
diagram of the physical model is shown in Fig. I. 

Heat Pipe Wall 

In this region, the progress can be modeled by the 
heat conduction equation in a hollow cylinder. The 
governing equation are expressed as follows: 
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Fig. 1 Schematic diagram and coordinate configuration of 
the heat pipe. 


Wick Stracture Region 

Initially, the working substance is in the soiid state. 
When heat is "added to the evaporator. ;he frozen working 
substance in the heated region is melted so that the liquid 
and solid states of the working substance exist in the wick. 
Fluid motion in the liquid region may then occur due to 
vaporization and condensation of the working substance. 
The liquid flow in the wick is considered to be unsteady 
two-dimensional incompressible laminar flow with 
negligible body forces. The fluid and wick structure are 
assumed to be in local equilibrium and the velocities in the 
axial and radial directions are the local area— averaged 
velocities over a cross section of a finite element of the wick 
region instead of the pore velocity or actual velocity. Also, 
the wick is assumed to be isotropic and homogeneous. The 
governing equations for the wick region are formulated by 
using the principles of the conservation of mass, momentum, 
and energy. The viscous dissipation terms in the energy 
equation are neglected. 

The continuity, momentum and energy equations are 
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region due to the extremely small density of the vapor. The 
vapor pressure is low and the temperature and pressure 
gradients are large in the axial direction, so the vapor 
velocity may reach the sonic velocity, and supersonic vapor 
flow and a shock front may occur in the condenser. Thus, 
the effects of compressibility, friction at the liquid— vapor 
interface and dissipation in the vapor should be considered 
in the mathematical model. The vapor flow may be 
considered to be axi symmetric, compressible, unsteady 
laminar flow and the governing equations for this flow are 
formulated with negligible body forces and heat sources as 
follows. In cylindrical coordinates, the continuity, 
momentum and energy equations are: 






(ii) 
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When the porosity approaches unity (no wick structure 
exists), the permeability approaches infinity. Therefore, 
equations (2-5) can be reduced to the N T avier-Stokes 
equation for unsteady two-dimensional incompressible 
iaminar flow. For the steady state, equations (2-5) also 
approach a special case given by Hong et al. (IS85). The 
expressions for the porosity and the effective thermal 
conductivity of the screen wicks are given by Chang (1987). 

The governing equations (2—5). however, are not 
always applicable. For example, during the second period 
the hauid state of the working substance exists in the wick 
but evaporation of the working substance a; the interface is 
negligible. Also, the liquid layer is so thin that the effect of 
natural convection in the liquid region is neglected. Thus, 
there is no fluid motion in the liquid region so that only 
equation ( 5 ) without the second and third terms is 
aoplicable and equations (2—4) are useful after the second 
period when there is liquid motion in the wick. 

In addition to equations (2—5), coupling conditions at 
the liquid— solid interface are needed: 
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Vapor flow dynamics 

Initially, the entire working substance is in the solid 
state so that the vapor space may be nearly evacuated. As 
the temperature at the interface increases, the vapor density 
also increases. Continuum flow in the vapor space is 
considered to be established when the mean free path, A, is 
substantially less than the mininum dimensions of the vapor 
flow passage (Holman, 1981), e.g., 

K n = ^<0.01 (9) 

♦ 

The transition temperature, T , of the vapor corresponding 
to the given dimension of the vapor space is expressed by 
using the kinetic theory of gases as follows (Jang, 1988): 



When the vapor temperature is greater than that calculated 
by equation ( 10 ), continuum flow is assumed to be 
established in the vapor space. 

WTien continuum flow is established in the vapor 
space during the start-up from the frozen state, complex 
flow phenomena are encountered in the continuum flow 
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The numerical and analytical solutions of the above 
equations under steady state conditions for annular and 
conventional heat pipes were given by Faghri (1986) and 
Faghri and Parvani (1988). Transient results are needed for 
the start— up from the frozen condition. 


Initial and Boundary Conditions 

Initially, a uniform temperature for the wall, wick, 
and vapor regions is assumed which is less than the melting 
temperature of the working substance. So, the entire 
working substance in the wick is in the solid state. The 
heat flux, convection and radiation boundary conditions are 
applicable at the outer surface of the evaporator and 
condenser of the heat pipe. Both ends of the heat pipe are 
assumed to be insulated. When the liquid motion exists in 
the wick structure, the no— slip condition is used for the 
velocities at the interface between the heat pipe wall and 
the wick structure. 

The boundary conditions at the liquid-vapor 
interface change during the start-up process. In the first 
and second periods, free molecular flow is prevalent in the 
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vapor space so that heat transfer through the vapor space is 
negligible. Therefore, the adiabatic condition, is employed 
at the liquid— vapor interface. When continuum fiow is 
established in the vapor space, vaporization and 
condensation of the working substance at the interface is 
considered by using interface energy and mass balances. 
The vapor temperature at the interface is assumed to be the 
saturation temperature corresponding to the vapor pressure 
at the interface. The Clausius— Clapeyron relationship is 
used to obtain the saturation temperature. The continuity 
of the axial velocity and shear stress is also used. The 
pressure difference at the liquid— vapor interface induced 
from the action of the surface tension is given by using the 
Laplace and Young equation. This condition couples the 
liquid and vapor momentum equations at the liquid— vapor 
interface. 

While part of the vapor space is in the continuum 
flow regime, free molecular flow also exists in the rest of the 
vapor space. Even though the heat transfer through the free 
molecular flow region may be negligible, the boundary' 
conditions at the border of the two regions are needed to 
solve the governing equations for the continuum flow region. 
Since a large temperature gradient exists in the continuum 
flow region in the condenser during start-up, most of the 
vapor may be condensed at the interface. Thus, the vapor 
penetration in the free molecular flow region may be 
minimal or penetration may occur in the immediate vicinity 
of the interface between the continuum flow region and the 
flee molecular flow region. Also, the temperature in the 
region of flee molecular flow remains unchanged except in 
the vicinity of the continuum flow region due to the near 
vacuum. Therefore, an imaginary plane, which is adiabatic 
and normal to the axial direction, is assumed to divide the 
two vapor flow regions at the point of the transition 
temperature. The dividing plane moves towards the cooled 
end of the heat pipe as the location of the transition 
temperature at the iiquid— vapor interface moves. 

The axisymmetric condition along the centerline of 
the heat pipe.is used in continuum flow region. The no-slip 
condition and the adiabatic condition at both ends of the 
heat pipe is also used for the velocities and temperature. 
The pressure and density at both ends. of the heat pipe are 
unknown but the vapor velocity near the ends of the heat 
pipe is so low that the axial pressure and density gradients 
are assumed to be zero. All of the boundary conditions for a 
mathematical model of the heat pipe start-up from the 
frozen state is summarized in Table L 

SIMPLIFICATION OF THE MODEL 

The mathematical model, equations (1—14), for the 
start— up behavior of the liquid metal heat pipe described in 
the previous section includes most of the physical 
phenomena which may occur in the heat pipe. Therefore, 
this model is very complex to solve numerically. The effect 
of some physical phenomena may be negligible, so a 
simplified model is derived to predict the performance of the 
heat pipe during the start-up period. For this purpose, 
assumptions are made based on the characteristics of the 
heat pipe and previous studies. 

The density of the liquid state of the working 
substance is much greater than that of the vapor state, so 
the velocity of the working substance in the wick structure 
is small. The thermal conductivity of the liquid metal is 
large and the thickness of the wick region is very thin. It is 
then assumed that the effect of the liquid flow in the wick 
structure is negligible and the wick structure is saturated by 
the working substance. Thus, the heat transport through the 
wick structure and working substance is by conduction only 
but the phase change of the working substance is considered. 
Under these assumptions, the same governing equation is 
also applicable to the heat pipe wall and the wick structure 
by using the proper properties for each region. The 


governing equation is given as follows: 
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The expression of (pc^)- for the wick structure region is 
given by equation (6). The thermal conductivity, K-, can 

be the thermal conductivity of the wall material for the wail 
region. When the effective thermal conductivity for the 
wick reeion is calculated by using the expression given by 
Chang ^ 1387), the thermal conductivity for the working 
substance is substituted by the solid conductivity, liquid 
conductivity, or the average value of the solid and liquid 
conductivity of the working substance corresponding to the 
solid, liquid, or mushy region, respectively. 

The vapor flow is also simplified further from the 
two— dimensional model to a one— dimensional model since a 
previous study (Jang et a!., 1989) shows that the 
one— dimensional transient compressible model described the 
vapor flow dynamics as well as the two-dimensional model 
for simulated heat pipe vapor flow. The transient 
compressible one-dimensional continuity, momentum and 
energy equation are written as follows: 
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Since the governing equations are simplified, the 
corresponding boundary conditions are also modified to 
match the governing equations. 


NUMERICAL PROCEDURES 

The governing equations are separately solved for 
each region. When the coupling is implemented at the 
interface, iterations are needed. To reduce the amount of 
computer time, non-iterative schemes are employed for 
each region. The well-known alternating direction implicit 
(ADI) method is used for the heat pipe wall and wick, and 
the phase change of the working substance during start-up 
is modeled by using the equivalent heat capacity method 
(Hsiao, 1985). This method approximates the rapid change 
of the heat capacity over the phase change temperature 
range, which is an artificially defined finite temperature 
range, AT, instead of using the Dirac function. In the 
numerical calculation, this property is evaluated based on 
the nodal temperatures. The implicit Beam-Warming 
method is used for the vapor flow dynamics. The detailed 
numerical method for the vapor flow in the heat pipe is 
described by Jang et al. (1989). 
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After continuum flow exisrs in the vapor space, 
equation (15) should be coupled wi;k equations (16-18) by 
using :he sane boundary conditions at the interface. The 
coupling of the governing equations for the vapor region to 
those for the wall and wick regions wpuid be’ achieved by 
using the heat flux and the saturation temperature at the 
interface. However, the heat flux at this interface and the 
saturation temperature are initially unknown, so that these 
boundary conditions should be assumed and iterations are 
needed for each time step until the coupling conditions are 
satisfied along the interface. The numerical procedure used 
for coupling is as follows: 

1. It is assumed that the liquid— vapor interface 
temperature is the initial temperature for the first 
time step. 

2. Solve for the temperatures in the wall and wick 
regions. 

3. Calculate the heat fluxes, q, at each node of the 
liquid— vapor interface by using the temperatures 
(Tj_, and Tj) in the wick region. 



4. These heat fluxes are used as the boundary 
conditions at the interface to solve for temperatures 
in the wall and wick regions. 

5. Use the same heat fluxes to solve the vapor 
temperature and pressure for the same period at the 
wall and wick regions. Obtain the saturation 
temperature, T s , by using the Clausius-Clapeyron 

relationship. 

6. Calculate the new heat fluxes, q', at the interface by 
using the saturation temperature, T in the wick. 



7. Compare the new heat fluxes, q', with the old heat 
fluxes, q, at each node of the interface. 

8. If the difference between the new heat flux and old 
heat flux is within an acceptable range, repeat steps 
4 to 7 for the next time step. 

9. If the difference between the new heat flux and the 
old heat flux is not within the acceptable range, 
assume new guessed heat fluxes, q M , by using the 
relaxation method and repeat steps 4 to 7 until the 
comparison of the results is acceptable. 

q" = q-i-a(q' -q) (21) 

10. Repeat steps 4 to 9 until the temperatures reach the 
steady state. 

When the coupling of the governing equations is 
attempted, some physical characteristics of the two regions 
are considered. Since the density of the vapor is much 
smaller than that of the liquid, the volumetric specific heat 
3 

pCp (— 8.8 J/m K) for the vapor is much smaller than that 

(= 1040.1 KJ/m 3 K) for the liquid. Therefore, a difference 
between the transient response times of the vapor region, 
and the wall and wick regions exists. The time step for the 
vapor space should then be much smaller than that for the 
wall and wick regions. The governing equations for ;he wall 
and wick regions are solved for one time step by using the 
heat flux assumed at the interface, and then the governing 
equations for the vapor space are solved by using the same 
heat flux at the interface and smaller time steps for the 
same period as the wall and wick regions. 


RESULTS AND DISCUSSION 

The governing equations for the wall and wick 
regions and the vapor flow axe separately solved, and then 
are coupled at the liquid-vapor interface. Therefore, the 
numerical methods and algorithms can be separately tested 
and compared to the available data. Since the working 
substance changes phase from the solid state to the liquid 
state during the start-up period from the frozen state, this 
effect should be incorporated into the numerical model. The 
solidification of.sodium in a square region is chosen to verify 
the numerical model and algorithm for governing equation 
f 15). The results are in agreement with available data 
(Rathjen and Jiji, 1971). The transient one-dimensional 
model for the vapor flow dynamics in the heat pipe has 
already been verified by Jang et al. (1989). The combined 
model is used to predict the performance of the liquid— metal 
heat pipe. 

Simulation of the heat pipe start-up in the initial stages 

Previous experimental results (Deverall et al., 1S70) 
show that successful start -up of the frozen heat pipe greatly 
depends on the boundary condition at the outer surface of 
the evaporator and condenser. For successful start— up from 
the frozen state, the heat input and output should melt the 
working substance in the condenser and allow sufficient 
liquid to return to the evaporator. All of the previous 
experimental results show the wall surface temperatures, so 
that even for successful start— up the status of the working 
substance with time is unknown. Thus, the governing 
equations which described the start— up behavior of a 
cylindrical heat pipe during the first and second periods are 
solved numerically to examine the effect of the boundary 
conditions and to recommend the optimum boundary 
condition. 

The physical model has evaporator, adiabatic, and 
condenser section lengths. olOrf, 0rf.aad 0*2 m, respectively. 
The radius of the vapor space and the inner and outer radii 
of the heat pipe wall are 0.00585, 0.008, and 0.01 m, 
respectively. The material for the heat pipe wall and wick 
structure is stainless steel (AISI 316). Sodium initially in 
the solid state at the ambient temperature is used as the 
working substance. The melting temperature of sodium is 
371 K and the transition temperature of 680 K is obtained 
by using equation (10) for the radius of the vapor space 
(0.90685m) by iteration. Three different boundary 
condition cases are chosen for the outer surface of the heat 

2 

pipe. A uniform input heat flux of 50 KW/m and radiative 
heat output are used on the evaporator for all three cases. 
An emissivity of 0.9 and a radiation reference temperature 
of 293 K are employed. The boundary conditions at the 
condenser and adiabatic sections are changed while that at 
the evaporator remains the same for all three cases. The 
liquid— vapor interface is assumed to be adiabatic due to the 
free molecular flow in the vapor space during initial two 
periods. 

For case 1, only radiation is used in the condenser to 
reject heat. Figure 2 shows the temperature distributions at 
the heat pipe outer wail surface and the liquid— vapor 
interface. As heat is added in the evaporator, the 
temperature in the evaporator increases and at 20 seconds 
the working substance in the evaporator is in the liquid 
state. However, the temperature in the condenser is not 
changed from the initial temperature and in the adiabatic 
section the temperature in the region adjacent to the 
evaporator increases due to axial conduction. Additional 
heat input in the evaporator increases the temperature 
above the transition temperature in the evaporator, but the 
temperature in the condenser is still not changed so that a 
large temperature gradient exists. In the adiabatic section, 
part of the working substance (sodium) is in the liquid 
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state. When the heat input is continued at the evaporator, 
vaporization occurs at the interface in the evaporator. 
However, the working substance in most of the adiabatic 
section and condenser is in the solid state. Therefore, the 
heat input in the evaporator should be small to prevent 
cry— Dili of the wick structure while the working substance 
in the solid state is melted. Even though successful 
start-up may be possible for this case, the start-up 
progresses very slowly. 



Fig. 2 Temperature distribution at the outer wall surface 
and liquid— vapor interface of the heat pipe wall with 
time for case 1. 
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For case 2, 10 kW/nT is added in the condenser in 
addition to the radiative boundary condition to assist in the 
start— up of the frozen heat pipe. Figure 3 shows the 
temperature distributions at the outer wall surface and 
liquid-vapor interface. Since a small amount of heat is 
added in the condenser, the temperature in the condenser is 
raised above the melting temperature. However the 
temperature in the adiabatic section is still below the 



Heat Pipe Axial Length, m 

Fig. 3 Temperature distribution at the outer wall surface 
and liquid— vapor interface of the heat pipe wall with 
time for case 2. 


melting temperature. Thus, liquid in the condenser cannot 
How to evaporator until the working substance in the 
adiabatic section liquenes. The temperature in the 
adiabatic section increases relatively faster than that for 
case 1 due to heat transfer at both ends of the adiabatic 
section. The start— up period may be shorter than that for 
case 1. 



Fig. 4 Temperature distribution at the outer wall surface 
and liquid— vapcr interface of the heat pipe wall with 
time for case 3. 

Fin ail jr, the adiabatic section is used as part of the 

condenser and 10 kW/m" of heat is input in the condenser 
section. Figure 4 shows the temperature distributions at the 
heat pipe surface and the liquid— vapor interface and Fig. 5 
shows the surface temperature at different times. Even 
though a large temperature gradient still exists along the 
axial direction, the working substance is completely, melted 
in the entire heat pipe in 80 seconds. When vaporization 
occurs in the evaporator, the working substance can flow 
from the condenser to the evaporator to prevent dry-out of 
the wick structure in the evaporator. Thus, a relatively 
large amount heat can be added at the evaporator without 
dry— out so that the start— up period is expected to be much 
less than those of cases 1 and 2. 

Transient Heat Pipe Operation 

To simulate the coupling of the governing equation 
fpr the wall and wick to that for the vapor flow, the same 
physical heat pipe model is used except that the adiabatic 
section is eliminated. Sodium is employed as the working 
substance. In order to concentrate on the coupling problem, 
it is assumed that continuum flow is established in the 
entire vapor space and the working substance is in the liquid 
state. To yield these conditions, a uniform initial 
temperature of 800 K, which is greater than the transition 
temperature (680 K), is used for the wall, wick, and vapor 
regions. The external surfaces of the heat pipe wall at the 
evaporator and condenser are .exposed to a uniform heat flux 
2 

of 50000 W/m and a convective boundary condition, 
respectively. A reference temperature of 300 K and a heat 

transfer coefficient of 100 W/m K are used for the 
convective boundary condition at the condenser section. 
The time step for the vapor space should be much smaller 
than that for the wall and wick regions, so a time step of 0.1 
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Fig. 5 Temperature distribution at the outer wall surface of 
the heat pipe wail with different times for case 3. 

second is employed for the wall and wick regions and time 

o 

step of 0.1 x 10 second is used for the vapor flow. 

Twenty nodes in the radial direction and 160 nodes in the 
axial direction are used at the wall and wick regions. Also, 
160 nodes are employed along the vapor space. A relaxation 
factor of a = 0.00003 is used to obtain the new guessed heat 
flux. 

Figure 6 shows the temperature distributions at the 
outer wall surface, liquid— vapor interface, and the 
saturation temperature for a time of 0.3 second. The 
temperature distribution at the outer wail surface is uniform 
within each section, and near the border between the 
evaporator and condenser the surface temperature abruptly 
changes corresponding to the boundary conditions at the 
surface. The interface temperature matches well with the 
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saturation temperature which is evaluated by using the 
Clausius— Clapeyron relationship with the vapo* pressure. 
The saturation temperature decreases gradually, figure 7 
shows the heat flux distribution at the interface. The new 
heat flux calculated is converged to the old heat flux. The 
maximum difference between the two heat fluxes is about 
10%. Even though the heat flux at the surface is relatively 
uniform, the heat flux at the interface is not uniform. Also, 
the heat flux at the interface is mnch less than that at the 
surface. This implies that most of the energy is used to raise 
the wall and wick temperature at this moment. Figure 8 
shows the vapor temperature, pressure, velocity, and density 
distributions. The variation of the vapor temperature, 



Fig. 7 Comparison of new heat flux with old heat flux at 
the interface for time of 0.3 second during transient 
continuum flow. 

pressure, and density is small. Also,, a Mach number of M = 
0.027 is obtained at the exit of the evaporator. 





Fig. 6 Temperature distributions at the outer wall surface 
and liquid— vapor interface, and saturation 
temperature for time of 0.3 second. 


Fig. 8 Axial variations of temperature, pressure, density 
and velocity for time of 0.3 second during transient 
continuum flow. 


CONCLUSIONS 

The start-up process of the frozen heat pipe is 
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described based on the experimental results. A complete 
mathematical model is developed to predict the start-up 
behavior of the heat pipe from the frozen condition. The 
simplified model is used to obtain the numerical results. The 
numerical results during the first and' second periods show 
that the heat flux distributions ' for the evaporator and 
condenser should be chosen to melt the working substance in 
the condenser prior to vaporization occurring in the 
evaporator. A small amount of heat input at the condenser 
helps start-up and a high heat rejection at the condenser 
during the start— up should be avoided. The coupling of the 
two-dimensional transient model for the wall and wick to 
the one-dimensional transient model for the vapor flow is 
achieved when continuum flow exists in the vapor space. 
During the transient operation, the heat flux distribution at 
the interface is quite different from that at the surface. 
Efforts to improve the present model will continue. 
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Table 1. Summary of the boundary conditions for a Mathematical Model 
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